Oscillations of hot, young neutron stars: 
Gravitational wave frequencies and damping times 
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We study how the frequencies and damping times of oscillations of a newly born, hot proto-neutron star de- 
pend on the physical quantities which characterize the star quasi-stationary evolution which follows the bounce. 
Stellar configurations are modeled using a microscopic equation of state obtained within the Brueckner-Hartree- 
Fock, nuclear many-body approach, extended to the finite- temperature regime. We discuss the mode frequency 
behaviour as function of the lepton composition, and of the entropy gradients which prevail in the interior of 
the star. We find that, in the very early stages, gravitational wave emission efficiently competes with neutrino 
processes in dissipating the star mechanical energy residual of the gravitational collapse. 
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I. INTRODUCTION 

The birth of a proto-neutron star (PNS) in a core collapse 
supernova is a very difficult phenomenon to model, since it 
requires not only accurate descriptions of the micro-physics 
of the collapsing matter, in particular of neutrino transport and 
related processes, but also of the violent dynamical processes 
occurring in the contracting-exploding star, which need to be 
treated in the framework of general relativity (see |U] for a 
recent review). 

The description of the subsequent PNS evolution is also 
challenging, because a PNS is a hot and rapidly evolving ob- 
ject. The physical processes which contribute to the star cool- 
ing and contraction, such as nuclear and weak interactions 
and energy and lepton number transport by neutrino diffusion, 
have to be included in dynamical simulations. Thus, most 
simulations of gravitational core collapse to a PNS end shortly 
after the core bounce and the launch of the supernova explo- 
sion - typically after a few hundreds of milliseconds - and 
only a few dynamical simulations extend to the first minute of 
the PNS lifelM. 

In this paper we are interested in this latest phase of the 
PNS life, when shock waves, neutrino winds, convection in- 
stabilities, and accretion flows are no longer dominant and the 
star cooling and contraction proceed on timescales of seconds, 
so that the evolution is quasi-stationary. In particular, we want 
to study how the frequencies and damping times of the PNS 
quasi-normal modes of oscillation depend on the star inter- 
nal structure. The main motivation is that the oscillations of 
a newly born PNS may be associated to gravitational wave 
signals with sizeable amplitudes, and with frequencies lower 
than those typical of mature neutron stars. This would favour 
their detection by the next generation of ground-based inter- 
ferometers LIGO/Virgo and their future version (ET) . 

The available dynamical simulations of the post-bounce 
evolution of a PNS indicate that, typically, the star goes 
through the following main steps. After the core bounce, a 
shock wave propagates through the outer PNS mantle, leaving 
behind a low-entropy core in which neutrinos are trapped, sur- 
rounded by a low-density, high-entropy envelope. The mantle 
accretes matter from the outer layers and rapidly contracts, 



losing energy due to electron captures and thermal neutrino 
emission. The supernova explosion lifts off the stellar enve- 
lope and, in a few tenths of seconds, due to extensive neutrino 
losses, the lepton pressure decreases and the envelope con- 
tracts. At this stage the PNS radius is about 20-30 km; the 
subsequent evolution can be described as a sequence of equi- 
librium configurations; this quasi-stationary evolution is the 
phase of interest for us. 

Simulations show that the diffusion of high-energy neutri- 
nos (of the order of a few hundred Me V) from the low-entropy 
core to the high-entropy envelope, from which they finally es- 
cape with energies of the order of a few tens of MeV, gen- 
erates a large amount of heat within the star, producing tem- 
peratures up to several tens of MeV; as a result, the core en- 
tropy approximately doubles, whereas the entropy of the en- 
velope decreases. In a few tens of seconds the PNS becomes 
lepton poor, but it is still hot. The net number of neutrinos 
in the interior is low, but thermally produced neutrino pairs 
of all flavors are abundant, and dominate the emission; the 
star cools down and entropy gradients are gradually smoothed 
out, while the average neutrino energy decreases, and neutri- 
nos mean free path increases; after approximately one minute 
it becomes comparable to the stellar radius, and the star be- 
comes neutrino transparent. By this time, the temperature has 
dropped to 1-5 MeV (« 1-5 x 10 10 K) and the star has radi- 
ated off almost all of its binding energy, becoming what we 
call a neutron star (NS). 

This brief summary of the first minute of the PNS life is de- 
liberately imprecise, because the details of the evolution de- 
pend on the assumptions on which the simulation is based. 
For instance, in |2Ha]< where the first minute after the core 
bounce is considered, the evolution is treated as a sequence 
of quasi-stationary states: the thermodynamical variables and 
the lepton fractions are determined by solving evolution equa- 
tions (for instance Boltzmann's equation to model neutrino 
transport), whereas at each time-step the stellar structure is 
found by solving the Tolman-Oppenheimer-Volkov equations. 
In |5J, instead, all quantities are determined through a time- 
evolution core-collapse code, which has been extended in or- 
der to describe the first ~ 20 s of the post-bounce processes. 
However, there are features which are common to different 
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studies; indeed, starting from an initial configuration charac- 
terized by a low-entropy core and a high-entropy envelope, 
due to neutrino processes the star goes smoothly through the 
following phases: 

• the entropy increases in the core while decreasing in the 
envelope; 

• entropy gradients gradually smooth out while the star is 
still hot; 

• the star progressively cools down, and the overall en- 
tropy decreases; 

• the evolution ends in the "cold", zero entropy, neutron 
star configuration. 

We remark that both in and in lH] the equation of 

state (EOS) of baryonic matter is a finite-temperature, field- 
theoretical model solved at the mean field level. 

In this paper, instead, we employ a microscopic EOS 
obtained within the Brueckner-Hartree-Fock (BHF) nuclear 
many-body approach extended to the finite temperature 
regime, as we shall discuss in Sec. [HI] We stress that our aim 
is not to model the PNS quasi-stationary evolution; we rather 
want to explore how the frequencies and damping times of the 
PNS quasi-normal modes depend on the physical quantities 
which characterize the quasi-stationary configurations, which 
are essentially the entropy profile (which will appear to be the 
most important in this respect) and the lepton composition. 

PNS are expected to be rapidly rotating; however in our 
study we neglect rotation, since we are primarily interested 
in the effects of the thermal and chemical evolution on the 
star oscillation frequencies, and in comparing the results with 
those of previous works which use different EOSs to model 
the PNS. 

The article is organized as follows. In Sec.[II]we briefly ex- 
plain how to compute the complex values of the quasi-normal 
mode frequencies using the relativistic theory of stellar per- 
turbations. In Sec.|nl]the derivation of the equation of state of 
hot nuclear matter used to model the PNS evolution is shortly 
illustrated. In Sec. [IV] we discuss how the different stages of 
a PNS quasi-stationary evolution are simulated by construct- 
ing stellar configurations with appropriate entropy and lepton 
fraction profiles. In Sec. [V] we compute and discuss the stel- 
lar parameters and the quasi-normal mode frequencies for the 
various configurations. Conclusions are drawn in Sec.lVIl 



II. THE QUASI-NORMAL MODES OF NEUTRON STARS 
A. Stellar perturbations 

In order to find frequencies and damping times of the 
quasi-normal modes (QNMs) of a star, we need to solve the 
equations describing non-radial perturbations of a (spherically 
symmetric) star in general relativity, which we briefly recall. 

The perturbed spacetime metric is expanded in tensor 
spherical harmonics, as (we use geometrized units, assuming 



c = G= 1) 

ds 2 = -e v (l + ZH^Y^e 1 ™} dt 2 
+e X (\-/Hi m Y hn e m, ^dr 2 
-2i(Or l+1 H[ m Y tm e imt dtdr 

+r 2 (l - r e K em Y em e im ^j (d$ 2 + sin 2 9 d<p 2 ) , (1) 

where CO is the frequency, F#„(0,<p) are the scalar spherical 
harmonics, and H^ m {r), K hn (r) describe the metric perturba- 
tions with polar parity, i.e., those transforming as (— 1) under 
a parity transformation. In this paper we do not consider per- 
turbations with axial parity, which transform as (— The 
functions y/(r), X(r) describe the unperturbed metric, and are 
found by solving the Tolman-Oppenheimer- Volkov equations. 
The four-velocity of the generic fluid element is 

«" = hS + Sh" = (e^ 2 , 0,0,0) +icoe-^ 2 (0, , 

(2) 

where ^ is the fluid element Lagrangian displacement, ex- 
panded in vector spherical harmonics 

Ut,r,#,<p) = -/V lm {r)dvY tm {V,<py W , 
^(t,r,^,q>) = -^(r)^ m (0,<p) e ifl » ■ (?) 

The fluid is also characterized by its pressure and energy 
density 

p(r) + Sp(t,r,-&,(p)=p(r)+/Sp em (r)Y em (V,<py' ot , 
e(r) + 8e(t,r,#,<p)= £ (r) + r e 8e em {r)Y bn ( , q> )e ia " . (4) 

We denote with 8 the Eulerian perturbations, and with A the 
Lagrangian perturbations, so that for instance the Lagrangian 
perturbation of the pressure is 

Ap = 8p + ^, (5) 
or 

i.e., 

Ap^ = 8p im + - W lm %. (6) 

r or 

Einstein's equations, linearized in the perturbations, yield a 
system of ordinary differential equations for the perturbed 
functions. Different equivalent sets of equations have been 
derived in the literature, using different gauge choices or dif- 
ferent combinations of the relevant equations ^- [T3l . In this 
pap er we use the formulation of Lindblom and Detweiler 
lIlOl [Till , consisting of a system of four first-order differen- 
tial equations (hereafter, the LD equations) for the functions 
{H( m ,K em ,W im ,X em }, where 

and algebraic relations which allow to compute the remain- 
ing functions {H^^H^ , V } in terms of the others (see Ap- 
pendix[S). To close the system, an EOS, relating the energy 
density e and the pressure p, has to be assigned. 
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B. The quasi-normal mode frequencies 

A QNM is a solution of the perturbation equations, which 
is regular at the center, continuous on the surface, and which 
behaves as a pure outgoing wave at infinity. Since in general 
relativity a non-radial oscillation is associated to gravitational 
wave emission, such solutions belong to complex frequencies: 



CO 



1 



(8) 



where a — 2nv and v is the pulsation frequency; Tqw is the 
damping time of the mode due to gravitational wave emis- 
sion. If the mode is unstable, its imaginary part is negative 
and — Tgw is the growth time of the instability. 

The procedure to find the QNM frequencies is the follow- 
ing: (i) We choose a value of / and a complex value of CO 
(since the background is spherical, the equations do not de- 
pend on the index m). (ii) We integrate Eqs. ( lAU . by impos- 
ing that the solution is regular at the center and that Ap = 
at the stellar surface [Eqs. 1A3MA4V 1. (iii) We impose that 
the solution and its first derivative are continuous on the stel- 
lar surface, and find the metric perturbations outside the star, 
(iv) In vacuum, the perturbed equations reduce to a simple, 
second-order differential equation (the Zerilli equation JA6b ), 
which we integrate up to radial infinity, (v) We check whether 
the solution satisfies the outgoing wave boundary condition at 
infinity ( |A8t which identifies a quasi-normal mode; we then 
repeat the procedure for different values of CO. The values of co 
which satisfy the outgoing wave condition can be found using 
a Newton-Raphson method. 

The polar QNMs are classified, following a scheme intro- 
duced by T.G. Cowling in Newtonian theory [ 13], on the basis 
of the restoring force which prevails when the generic fluid 
element is displaced from the equilibrium position. Thus, 
we have a g-mode if the restoring force is mainly provided 
by buoyancy or a p-mode if it is due to a pressure gradient. 
The frequencies of the g-modes are lower than those of the p- 
modes, and the two sets are separated by the frequency of the 
fundamental (/-) mode, which is related to a global oscillation 
of the fluid. In general relativity there exist further modes, 
named w-modes [14], that are purely gravitational, since they 
barely excite the fluid motion. Other classes of modes are as- 
sociated to NS features which are not included in the present 
model, like rotation, magnetic fields, the crust rigidity. 



C. Sound speed 

A neutron star at the end of its evolution is cold and isen- 
tropic, matter is in beta-equilibrium and can be described by a 
barotropic EOS p = p(s). Conversely, when the star is young 
and hot the EOS cannot be expressed in a barotropic form, 
since the pressure depends non-trivially on the entropy and on 
the composition, i.e., 



p = p(s,s,Xi) 



(9) 



Therefore, to solve the perturbed equations the profiles of en- 
tropy and particle fractions, respectively s(n) and X/(n), are 



also needed. In Eq. |9) n is the baryon number density, 
s=S/A is the entropy per baryon, and X\ = fii/n is the fraction 
of the j'-th particle. Usually matter is locally in beta equilib- 
rium and neutrinos are trapped, therefore the dependence on 
the composition {x,} reduces to a dependence on the lepton 
fraction Y e = x e + x Ve only. 

The perturbed equations ( IA1 b depend explicitly on the 
sound speed c%, which relates the Lagrangian perturbations 
of pressure and energy density, 



A P : 



: c?Ae . 



(10) 



c s is defined as the following thermodynamical derivative 



^ e J adiabatic 



(ID 



where "adiabatic" means that the derivative is performed 
keeping fixed the entropy and the fractions of those particle 
species which do not change during the pulsation [15]. 

To clarify this statement, let us consider a fluid element os- 
cillating with period f osc about the equilibrium position. The 
following equation holds: 



Ap- 



dp 
ds 



Ae 



dp 
ds 



As- 

E,Xi 



dp 

dxi 



Axi 



(12) 



Since we are considering adiabatic perturbations, the fluid ele- 
ment does not exchange heat with its surroundings and As = 0. 
Furthermore, the displaced fluid element has a composition 
different from the surrounding fluid even though nuclear reac- 
tions, acting on a timescale f re act> tend to eliminate this differ- 
ence. The two limiting cases are: 

i) Wet ^ fosc; in this case the fluid element composition 
does not change during the oscillation, i.e., Ax, = 0, and by 
combining Eqs. ( fTOb and dT2b we find 



dp 
de 



(13) 



ii) freact "C t osc ; the fluid element composition changes, be- 
coming that of the surrounding fluid. By replacing the com- 
position profile x, = Xi(e,s,Y e ) in ©, it is possible to express 
the EOS as p = p(e,s,Y e ). Eqs. ( TTOt and (fT2l then give 



dp 
de 



(14) 



s,Y e 



for beta-stable, neutrino-trapped matter. 

For the PNSs in quasi-stationary evolution we consider in 
this paper, typical oscillation periods are of the order f osc ~ 
10~ 3 s, while weak interactions timescales are lfl6ll 



5 x 10 6 s 
(T/10 9 K) 6 



20 s 



(T/1Q 9 K) 4 



(15) 



for modified and direct Urea processes, respectively. In the 
first seconds of a PNS life T « (1 -4) x 10 11 K, thus f^t < 
fosc- Therefore the stellar pulsations always occur in local beta 
equilibrium. Timescales of strong nuclear reactions are even 
smaller. 
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D. Some considerations on the /- and #-modes 

For old, cold neutron stars, the frequency of the /-mode, 
Vf, is in the range 1-3 kHz, and the damping time, TV, is of the 
order of a few tenths of seconds. According to the Newtonian 
theory of stellar pulsations, V/ scales as the square root of 
the star average density, and this behaviour is maintained in 
the relativistic theory, according to which the damping time 
scales as x f ~ R 4 /M 3 CMH. 

The g-modes are directly related to the thermodynamical 
properties of the star. Indeed, their presence can be traced 
back to the Schwarzschild discriminant 120], 



(16) 



dr dr dr \ cin 



where c 2 Q = d fJ^ r . The radial acceleration of a fluid element 
displaced from equilibrium by Ar is 



-A/2 



(e+p) 



2 n 2 



dp 
~dr~ 



S(r)Ar. 



(17) 



Therefore, if S(r) > the fluid element oscillates about the 
equilibrium position, whereas if S(r) < it is accelerated 
away from equilibrium. It follows that, if in some region of 
the star S(r) < 0, this region is convectively unstable and the 
star admits a set of unstable g-modes, otherwise the g-modes 
are stable. If S(r) vanishes identically, the star does not have 
g-modes (all g-modes degenerate to zero frequency). This is 
the case if the neutron star is cold and old, since the EOS 
is barotropic and c 2 = Similar information is contained 
in the Brunt-Vaisala frequency which, in a relativistic frame- 
work, is defined as I2II1 



N 2 (r) = 



¥. 



(e + P )c 2 2 



(18) 



It has been shown that, although the Brunt-Vaisala frequency 
changes by many orders of magnitude throughout the star, it 
allows to estimate some g-mode frequencies of Newtonian 
stars. For instance, in 12211 the frequency of higher-order g- 
modes of main sequence stars is computed using the following 
formula 



(lK + l + n e 



R dr \m 



(19) 



where R is the stellar radius, K is the order of the g-mode, and 
n e is the effective polytropic index of the outer layers of the 
star. However, the Brunt-Vaisala frequency cannot be used to 
estimate neutron star g-mode frequencies 112 111 : these frequen- 
cies can only be found by solving the perturbation equations, 
as we do in the present paper. Nevertheless the following con- 
siderations will be useful to interpret the results we will show 
in the following. Eq. ( fT9] l indicates that higher frequency val- 
ues correspond to larger values of S(r) (i.e., of |iV(r)|) inside 
the star. Since we shall assume dY e /dr = (see Sec. HVl i. we 
have 



S(r) = 



dp 



dr 



dp 
de 



s,Y t 



de 
dr 



dp 
ds 



ds 
dr 



(20) 



we may expect that higher g-mode frequencies correspond to 
larger entropy gradients. 



E. The damping time of quasi-normal modes 

A QNM of a PNS is characterized by the pulsation fre- 
quency and by the damping time Tgw- Its value is important 
because it shows how fast the pulsation energy can be dis- 
sipated through gravitational wave emission, and it must be 
compared to the timescale Tdiss associated with other dissipa- 
tive processes which may compete with GW emission. These 
include viscosity, heat transport, neutrino diffusion, etc. (phe- 
nomena which we are neglecting in our model). In the first 
minute of life of a PNS, T diss ~ 10 - 20 s JHIIl (see also 
the discussion in |0]). Thus, if Tqw T<iiss> the oscillations 
are mainly damped by gravitational wave emission, and vicev- 
ersa. We also remark that, if a QNM is unstable, the instability 
can grow only if Tgw <C Tdi ss - 

As long as Tgw <C Tdi ss , when a star oscillates in a QNM, 
the pulsation energy changes in time as J9)] 



£puls(0 



J puls 



(o), 



,-2r/T GW 



and the power radiated in gravitational waves is 



^puls 



2£ P uis/ Tgw ■ 



(21) 



(22) 



Thus, smaller QNM damping times are associated with a more 
efficient gravitational wave emission. In the case of cold NSs, 
the /-mode has always the smallest damping time, but this is 
not always the case for PNSs, as we shall discuss later. 

Although we shall compute the damping times of all modes 
by direct integration of the perturbed equations, it is useful to 
give some approximate formula which will allow us to explain 
some results of the next sections. From Eq. (l22l we find 



tqw ~ 2£puis/£GW • 



(23) 



The (approximate) expressions of £p U i s and Lqw (the latter 
is obtained using the quadrupole formalism) can be found in 
j2lll26Tl and are, in terms of the perturbation functions defined 
in this paper, 



£puls ~ ^ ° L J Q dr r ( £ + 

x \\W em \ 2 + £(£+l)\V 



and 



4n , 

£gw ~ 



dr r de 



.4 2r.tm 



where 



8e lm = -r e 



-V/2 



e- x > 2 
X im + e r W 



fall 



(24) 



(25) 



(26) 
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HI. THE EQUATION OF STATE OF HOT NUCLEAR 
MATTER 



A. BHF many-body approach 

One of the most advanced microscopic approaches to the 
EOS of nuclear matter is the Brueckner theory B27I1 . recently 
extended to the finite-temperature regime within the Bloch- 
De Dominicis formalism [28]. In this approach, the essen- 
tial ingredient is the two-body scattering matrix K, which, 
along with the single-particle potential U, satisfies the self- 
consistent equations 



(12|#(W)|34) = (12|V|34) 
+ Re£(12|V|3'40 [1 - nF(3 ' )][1 

3', 4' 



7 (4')] 



W-E* 



l£ 



{3'4'\K(W)\34) 



and 



U(l)=£n F (2)(l2\K(W)\12) A , 



(27) 



(28) 



where 1,2, ... generally denote momentum, spin, and isospin. 
Here V is the two-body interaction, W = E\ + Eo repre- 
sents the starting energy, and Ei = fc?/2m; + U(kj) the single- 
particle energy; n (k) is the Fermi distribution at finite tem- 
perature. For assigned partial densities and temperature, 
Eqs. dZTt and d28l have to be solved self-consistently along 
with the following equations for the auxiliary chemical poten- 
tials |Uj, 



i 



,P(E,(k)-ji,) + i 



(29) 



and the baryon number density is n = n i- 

At finite temperature the EOS, and all thermodynamical 
quantities, can be obtained from the grand-canonical poten- 
tial density CO. In the Bloch-De Dominicis approach, (0 can 
be written as the sum of a mean-field term and a correlation 
contribution 127, 291, 



0) : 



E 

k 

1 



1 



In (l+e 



-^+n F (k)U(k) 



Z,n F (k)U(k). 



In this framework, the free energy density is 

f = co + nfi, 



(30) 



(31) 



and all remaining thermodynamical quantities of interest, 
namely, the "true" chemical potential ji, pressure p, entropy 
per baryon s, and energy density e can be computed from it as 



an 



p = n 



)n 

2 d(f/n) 



dn 

= -121 

S ndT' 

£ = f + Tns + m„n 



= lM-f, 



(32) 

(33) 

(34) 
(35) 



(m„ neutron mass). Since at zero temperature the non- 
relativistic microscopic approaches do not correctly repro- 
duce the nuclear matter saturation point, no ~ 0.17 fm -3 , 
E/A w —16 MeV, three-body forces (TBF) among nucle- 
ons are usually introduced in order to correct this deficiency. 
Given the current lack of a complete microscopic theory of 
TBF, we have adopted the phenomenological Urbana model 
13011 . which consists of an attractive term due to two-pion ex- 
change with excitation of an intermediate A resonance and a 
repulsive phenomenological central term. For simplicity, we 
reduce the TBF to a density-dependent two-body force by av- 
eraging over the position of the third particle, assuming that 
the probability of having two particles at a given distance is 
reduced according to the two-body correlation function. The 
corresponding EOS at zero temperature reproduces the nu- 
clear matter saturation point correctly i3l1 - [33ll . and fulfills 
several requirements from the nuclear phenomenology. In all 
calculations presented in this paper we use the Argonne Vig 
nucleon-nucleon potential [34] together with the phenomeno- 
logical Urbana TBF. 

Results for symmetric nuclear matter and purely neutron 
matter have been obtained for different values of the temper- 
ature, and are discussed in Il35l - [38ll . In particular, in Ref. 113711 
useful numerical parametrizations of the EOS are given that 
are employed in the current work. 

The Brueckner approach provides a realistic modeling of 
nuclear matter only at densities above about half normal nu- 
clear matter density. Below this threshold, clusterization sets 
in, and the system becomes inhomogeneous. Therefore, in 
this "low-density" regime another theoretical approach has to 
be used, and we employ the EOS of Shen IB9I1 . which is es- 
sentially a liquid-drop-type model at finite temperature. 

Of course, since two different theoretical descriptions of 
the same state of matter are involved, the joining of the two 
EOSs requires the thermodynamical observables to be contin- 
uous functions of the baryon density. In practice we perform a 
Maxwell construction by equating pressure and chemical po- 
tentials of the low- and high-density sectors, and verify that 
the other thermodynamic variables do not exhibit significant 
discontinuities at the transition point. In this way, a very wide 
range of baryon density is spanned, from the iron density at 
the surface up to 8-10 times the nuclear saturation density in 
the core. Further details are discussed in the following sub- 
section. 



B. Composition and EOS of hot stellar matter 

In neutrino-trapped j3 -stable nuclear matter, the chemical 
potential of any particle i = n,p,l is uniquely determined by 
the conserved quantities, baryon number B,-, electric charge 
Qi, and weak charges (lepton numbers) L. , L^: 



(36) 



For stellar matter containing nucleons and leptons as relevant 
degrees of freedom, the chemical equilibrium conditions read 
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explicitly as 



entropy profile 



ju„ -ii p = pL e - ju Ve = Mm + Mv M 



(37) 



At given baryon number density n, these equations have to be 
solved together with the charge neutrality condition 

£fi« = 0, (38) 
and those expressing conservation of lepton numbers 



Y,=xi -Xf 



l V; 7 



(39) 



When neutrinos have left the system, their partial densities and 
chemical potentials vanish and the above equations simplify 
accordingly. We fix the muon fractions to = 0, and let Y e 
assume a finite value different from zero in neutrino-trapped 
matter. 

The nucleon chemical potentials are obtained from the free 
energy density /, Eq. d3TT i. 



ant 



(40) 



and the chemical potentials of the non-interacting leptons are 
obtained by solving numerically the free Fermi gas model at 
finite temperature. Once the hadronic and leptonic chemical 
potentials are known, one can proceed to calculate the compo- 
sition of the /3 -stable stellar matter, and then the total pressure 
p through the usual thermodynamic al relation 



p = n 



,d {f/n) 
dn 



(41) 



An important feature of the low-density domain is the 
treatment of neutrino trapping. Physically, neutrinos escape 
rapidly from the low-density matter during the PNS evolution, 
and the lepton number is not conserved anymore. This effect 
can be roughly modeled by introducing a neutrino-sphere in- 
side which neutrinos are trapped. Typical model-dependent 
values for the location of the neutrino-sphere found in the 
literature are 2 x 10" 3 fin -3 H, 6 x 1CT 4 fm~ 3 EH, and 
2 x 10~ 5 fm~ 3 14211 . Given these variations, we choose the 
following "natural cutoff" procedure: when imposing a con- 
stant Y e at any density, at a certain threshold number den- 
sity n v ps 10~ 5 — 10~ 6 fm~ 3 , the electron fraction x e becomes 
equal to Y e , and neutrinos disappear naturally. For lower den- 
sities we consider the matter untrapped. This simple proce- 
dure avoids making assumptions about the neutrino-sphere, 
although a more satisfactory treatment of neutrino trapping is 
required; but this is beyond the main goal of the present paper. 



IV. PROTO-NEUTRON STAR STELLAR MODELS : 
ENTROPY AND LEPTON FRACTION PROFILES 

We shall now construct equilibrium stellar models, all with 
a fixed baryonic mass Mg = 1 .5Mq (a conserved quantity dur- 
ing the stellar evolution), and with different entropy and com- 
position profiles. These configurations will be used to simu- 
late the quasi-stationary evolution of a PNS, and to compute 
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FIG. 1. (Color online) The profiles of entropy (upper panel) and 
temperature (lower panel) are plotted versus the enclosed baryonic 
mass for the models -Pi ,5,0.35 (ft = 30.3 km) and £^,3,0.30 (ft — 16-5 
km). 



how the stellar parameters and the quasi-normal mode fre- 
quencies change during the evolution. As discussed in the in- 
troduction, the quasi-stationary evolution typically starts with 
configurations characterized by a low entropy per baryon in 
the core (order of s ~ 1 at the center, see for instance J2(]) 
and a large entropy per baryon in the envelope (order of s ~ 5 
or larger). Thus, we shall consider as "initial" the configura- 
tion with an entropy per baryon profile made of two constant 
pieces, s c = 1 in the core and s e = 5 in the envelope, with a 
smooth junction between them. Furthermore, as discussed in 
Sec. IIII Bl to model neutrino trapping we shall assume that 
the lepton fraction Y e is constant throughout the star (up to a 
threshold density, below which Y e — x e ). 

As long as the evolution proceeds, entropy gradients are 
gradually smoothed out: the core entropy increases, the enve- 
lope entropy decreases, neutrinos escape from the surface and 
the star progressively cools down. To model this evolution, we 
construct EOSs and stellar configurations corresponding to in- 
creasing values of s c , decreasing values of s e and decreasing 
lepton fraction; then, to decreasing values of both s c , s e , and 
decreasing lepton fraction. Each configuration depends on the 
three constants s c ,s e ,Y e , and it is labeled by P Sc>Set Y e . 
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TABLE I. Stellar models with fixed baryonic mass Mg = 1.5Mq 
corresponding to different entropy profiles and lepton fractions. The 
gravitational mass M, radius R, central temperature T c , and central 
neutrino fraction x c v are tabulated for each profile. 



Sc 


s e 


Ye 


M/Mq 


R(km) 


T c (MeV) 


A v 


1.0 


5.0 


0.38 


1.43 


31.5 


19.8 


0.052 


1.0 


5.0 


0.35 


1.43 


30.3 


20.2 


0.041 


1.0 


5.0 


0.28 


1.42 


29.4 


21.2 


0.020 


1.5 


4.5 


0.33 


1.43 


24.6 


30.3 


0.035 


1.5 


4.5 


0.32 


1.42 


24.4 


30.5 


0.031 


2.0 


4.0 


0.32 


1.43 


21.5 


40.5 


0.033 


2.0 


4.0 


0.30 


1.43 


21.3 


41.0 


0.027 


2.0 


3.0 


0.30 


1.42 


16.5 


41.2 


0.026 


2.0 


3.0 


0.28 


1.41 


16.4 


41.8 


0.021 


2.0 


2.0 


0.28 


1.41 


14.5 


41.6 


0.020 


2.0 


2.0 


0.23 


1.40 


14.1 


42.9 


0.010 


1.0 


1.0 


0.23 


1.37 


12.5 


20.2 


0.007 


1 


1 


x v =0 


1.36 


12.2 


20.9 


0.000 




T = 


= 


1.35 


11.9 


0.00 


0.000 



To describe the latest stages of the PNS evolution we also 
consider two constant entropy profiles: (i) one with s c = s e = 
1 and no neutrino trapping, Pi.i. Xv= o, with Y e varying from 
0.10 at the center to 0.44 at the stellar surface; (ii) a zero- 
temperature profile with no neutrino trapping, Pt=q, which 
describes a cold, old NS, with electron fraction varying from 
0.09 at the center to 0.44 at the surface. 

We show in Fig. [T]the profiles of the entropy per baryon 
(upper panel) and of the temperature (lower panel) as a func- 
tion of the enclosed baryonic mass, for the models ^1,5,0.35 
(radius R = 30.3 km) and ft, 3,0.30 (radius R = 16.5 km). In 
order to avoid sharp transitions from the core to the envelope 
region, we adopt a cubic interpolation for the entropy between 
the two regions. Thus, the entropy is a continuous function of 
the density. However, as a consequence of the Maxwell con- 
struction used to join the Shen EOS (low-density region) to 
the BHF EOS (high-density region), as discussed in the previ- 
ous section, there is a weak discontinuity both in the entropy 
and in the temperature profile, when plotted as a function of 
the enclosed mass as in Fig. [TJ We have checked that the re- 
sults presented in the next section are not influenced by these 
discontinuities. 

In Table U we show the quantities which characterize the 
stellar models associated to different profiles, namely gravita- 
tional mass, stellar radius, temperature and neutrino fraction 
at the center of the star. The dependence of the stellar param- 
eters on the temperature and lepton fraction profiles will be 
discussed in the next Section. 



V. RESULTS 

In this section we discuss the behaviour of the stellar radius 
and of frequencies and damping times of the QNMs com- 
puted for stellar models with different entropy profiles and 
lepton/neutrino fraction content, in order to understand how 
these quantities are affected by the PNS internal structure. 



TABLE II. Frequencies (in Hz) and damping times (in s) of the 
QNMs g\,f,pi for different stellar models with baryonic mass 
1 .5 Mq , lepton fraction Y e = 0.32, and different entropy profiles. The 
central temperature T c (in MeV) and the stellar radius R (in km) are 
also shown. 





s e 


As 




R 






v f 


V 


v ,n 


V, 


1.0 


5.0 


4 


20.6 


29.6 


906 


6.27 


1194 


4.42 


1528 


0.75 


1.5 


4.5 


3 


30.5 


24.3 


910 


42.9 


1346 


0.76 


1845 


0.55 


1.0 


4.0 


3 


20.2 


18.4 


870 


793 


1741 


0.27 


2574 


0.99 


2.0 


4.0 


2 


40.5 


21.5 


669 


2xl0 3 


1449 


0.45 


2097 


0.72 


2.0 


3.0 


1 


40.7 


16.8 


492 


6 x10 s 


1714 


0.25 


2977 


1.64 



Let us consider the dependence on the entropy profile first, 
and fix the value of the lepton fraction to Y e = 0.32. We 
compute and compare the mode frequencies and damping 
times of the following stellar configurations: Pi 5 o.32> Pi,4,0.327 
Pi.5,4.5,0.32> P2,4,0.32> an d P2,3,o.32- The core-envelope 'entropy 
jumps' are As = 4,3,3,2, 1, respectively. In Table Ull we show 
for each profile the central temperature, the radius, and the 
frequencies and damping times of the QNMs gi,f,pi- These 
data allow us to discuss how the different quantities change 
with the entropy profile. 

As a general rule, the radius is larger if the star is hotter, 
or equivalently, if it has a larger entropy per baryon. This is 
indeed confirmed comparing for example the profiles Pi 5 0.32 
and Pi 4 o.32- The temperature (entropy) at the center is the 
same, but the first model has larger entropy and temperature 
in the envelope; its radius, R = 29.6 km, is larger than that of 
the second model, R — 18.4 km. This behaviour is confirmed 
by comparing P2,4,o.32 an d P2,3,0.32- In a similar way, the radius 
depends on the entropy in the core, although the dependence 
is weaker, because the envelope has more freedom to expand 
than the core; for instance, the configuration P24 o.32 nas a 
radius R — 21 .5 km, larger than R = 18.4 km of Pi ,4,0.32- 

As discussed in Section III Dl the frequency of the first g- 
mode depends mainly on the core-envelope entropy jump: 
higher values of As = s e — s c correspond to larger g-mode fre- 
quencies. Furthermore, as argued in lEBTl . the g-mode fre- 
quency has also a (weaker) dependence on the central temper- 
ature; indeed, the configurations P1.5 4.5 0.32 an d Pi 4 0.32 nave 
the same entropy jump, but the former has a larger central 
temperature T c and larger g-mode frequency. 

Table [TT] shows that, as the entropy jump decreases, the 
damping time of the first g-mode increases dramatically: for 
Pi,5,0.32 it is T gi = 6.27 s, for P2,3,o.32 it is T gl = 6x 10 5 s. 
Since the evaluation of the damping time is more sensitive to 
the numerical procedure than that of the mode frequency, we 
have also computed this quantity by using the approximate 
formula given in Eq. d23l ), and indicate its value as r| st ; we 
find a reasonable agreement with the data of Table [TT1 

As=4: P pu i s = 1.047km, L GW = 0.377, ^ st = 5.55s 
As=l: P puls =0.111km, L GW = 3.8 x 10" 6 , Tf l = 6xl0 4 s. 

Thus, the sharp increase of T gl as As decreases is due to 
the sharp decrease of the gravitational luminosity Low [see 
Eq. d23H . As shown in Eq. d25t . Low is (modulo a numerical 
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FIG. 2. (Color online) Comparison of the perturbed energy density 
8e(r), Eq. J26l >. (upper panel) and of the function /gwM given in 
Eq. J42b (lower panel) for the g-mode of the stellar configurations 
^1 5,0.32 an d ^2. 3,0.32- Both functions have been normalized in such a 
way that the mode pulsation energy is E vu \ s = 1 km. 

factor) the squared integral of the function 

/gw - vV5e*" , (42) 

where v is the mode frequency and 8e is the perturbed energy 
density for the considered mode. 

In order to understand why Low decreases so much when 
the entropy jump decreases, we plot in Fig.|2]5e (upper panel) 
and /gw (lower panel) as functions of r, for the stars with pro- 
files Pi 50.32 an d Pi, 3,0. 32- It is obvious that due to the larger 
radius of the former configuration and the presence of the fac- 
tor r 4 in Eq. d25l ), the emitted power Lqw is much larger (in- 
deed the main contribution comes from the envelope), and the 
damping time is strongly reduced. In addition, the g-mode fre- 
quency is also larger in the former configuration, and this con- 
tributes further to a larger gravitational wave emission, since 
7 G w ~ V 3 . 

Let us now consider the /-modes. For a cold neutron star, 
the /-mode frequency scales as the average density of the star 
and the damping time scales as iy ~ /? 4 /M 3 . From the data of 
Table HT1 we see that V/ increases as the radius decreases (the 
gravitational mass is nearly the same for all configurations), 



TABLE III. Frequencies (in Hz) and damping times (in s) of the 
QNMs gi,f,pi for stellar models with baryonic mass 1.5 Mq, en- 
tropy per baryon in the core s c = 1 and in the envelope s e = 5, and 
different values of the lepton fraction Y e . The radius of the star (in 
km) and its gravitational mass (in solar masses Mq) are also shown. 



Y e 


R 


M 


v« 




v f 








0.38 


31.5 


1.43 


863 


6.78 


1116 


9.75 


1415 


1.00 


0.36 


30.6 


1.43 


883 


6.62 


1147 


6.83 


1463 


0.89 


0.32 


29.6 


1.42 


906 


6.25 


1194 


4.44 


1527 


0.75 


0.30 


29.4 


1.42 


910 


5.99 


1209 


4.01 


1543 


0.73 


0.28 


29.4 


1.42 


908 


5.71 


1216 


3.96 


1546 


0.72 



while the damping time decreases; however, both V/ and xy 
do not follow quantitatively the cold star scaling laws. The 
first p-mode frequency has a behaviour similar to that of the /- 
mode, whereas the damping time seems to be quite insensitive 
to changes of the entropy profile. 

Finally, we consider a sequence of stellar models with a 
fixed entropy profile, i.e., s c = 1 in the core and s e = 5 in 
the envelope, and lepton fraction varying in the range Y e = 
0.38, . . . ,0.28. The frequency and the damping times of the 
gi, /, and pi-modes are shown in Table Hill together with the 
radius and the gravitational mass of the star. From these data 
we see that the star radius is a slightly decreasing function of 
the lepton fraction, and that the behaviour of the /- and p\- 
frequency as a function of the star radius is similar to that 
described above. Overall, the data show that the dependence 
of the QNMs eigenfrequencies on the lepton fraction is much 
weaker than that on the entropy profile. 

A. QNM eigenfrequencies and PNS quasi-stationary evolution 

As mentioned in the Introduction, numerical simulations 
show that in the early phases of a PNS life the entropy profile 
has a characteristic evolution which mainly depends on neu- 
trino diffusion processes, and which can be divided in three 
essential steps: 

1 . the entropy per baryon is initially (a few tenths of sec- 
onds after bounce) larger in the envelope and lower in 
the core; 

2. the entropy increases in the core and decreases in the 
envelope, reaching a roughly uniform profile; 

3. the entropy decreases throughout the star, which even- 
tually becomes a cold neutron star. 

The entire process takes about a minute, but we cannot as- 
sign precise temporal labels to each step, because they depend 
on the details of the initial conditions after the bounce and on 
the dynamical modeling of the evolution, which is beyond the 
scope of our work (an example of this evolution is shown in 
Fig.9ofRef.il). 

In this section we construct a sequence of stellar configura- 
tions, listed in Table [IV] which captures the main qualitative 
features of a PNS evolution described by steps 1 to 3. Each 
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TABLE IV. Frequencies (in Hz) and damping times (in s) of the 
QNMs gi,/,Pi for a sequence of stellar models which mimic the 
quasi-stationary evolution of a PNS with constant baryonic mass 
1 .5 M,-.)- The star radius (in km) is shown in column 5. 



i 


s c 


s e 


Y e 


R 






v f 








1 


1.0 


5.0 


0.35 


30.3 


890 


6.54 


1162 


5.89 


1484 


0.84 


2 


1.5 


4.5 


0.32 


24.4 


910 


42.9 


1346 


0.76 


1845 


0.55 


3 


2.0 


4.0 


0.30 


21.3 


667 


2.3 xlO 3 


1452 


0.44 


2125 


0.73 


4 


2.0 


3.0 


0.28 


16.4 


485 


7.6 xlO 4 


1717 


0.25 


3133 


1.80 


5 


2.0 


2.0 


0.23 


14.1 







1790 


0.23 


4134 


2.59 


6 


1.0 


1.0 


Xy = 


12.2 







1896 


0.21 


5879 


2.98 


7 




T = 


= 


11.9 







1898 


0.21 


6006 


3.52 



profile is labeled by a number i, which gives the ordering in 
time of the simulated evolution. Configurations from i = 1 to 
i — 4 (envelope entropy larger than core entropy) refer to the 
transition from step 1 to step 2, which ends with configuration 
5, for which the entropy distribution becomes uniform, but the 
star is still hot. Then it cools down (configurations 5 to 6) and 
ends as a zero-temperature NS reaching configuration 7 (step 
3). During this "evolution" the lepton number decreases. We 
have also considered a different sequence, in which the lepton 
fraction decreases "more rapidly," but the results are very sim- 
ilar to those obtained with the sequence shown in Table |IV| 

For each configuration we compute the frequencies and 
damping times of the QNMs gi,f,p\- Their values are given 
in Table [IVj and are plotted in Fig. [3] versus the number i 
which identifies the configuration as explained above. We 
remark that, as shown in Fig. [3] the gravitational damping 
time of the gi-mode sharply increases for i > 2, while the 
mode frequency sharply decreases. However, as discussed in 
Sec.HU as soon as z g becomes comparable to Tdiss ~ 10 — 20 
s, the mode becomes ineffective with respect to gravitational 
wave emission, since the stellar oscillations are damped by 
non-gravitational dissipative processes. 

The most interesting result which emerges from Table ITVl 
and Fig.[3]is that at earlier times, i.e., for i < 2, the frequen- 
cies of the gi,f,pi modes cluster in a small region around 
1 kHz, and then tend quite rapidly to the values appropriate 
for a cold NS (remember that our entire sequence should cover 
approximately a minute of the PNS evolution). This behaviour 
is similar to that found in where the quasi-stationary 
evolution sequence was obtained using a finite-temperature 
EOS derived within the mean field approach, treating neutrino 
transport using the diffusion approximation. The fact that in 
the very early stages the mode frequency is of the order of 
1 kHz (or lower) is important for gravitational wave detec- 
tion, because the sensitivity of ground-based interferometers 
LIGO/VIRGO decreases quite significantly at larger frequen- 
cies. 

Another interesting point to note is that during the early 
evolution the damping time of all modes is smaller than 10 s. 
This means that gravitational wave emission is effectively 
competing in removing energy from the star with dissipation 
processes related to neutrino viscosity, diffusivity, and thermal 
conductivity, since typical neutrino timescales are of the order 
of 10-20 s (see also section 2.1 of Q6[] for a detailed discussion 
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FIG. 3. (Color online) Frequencies (upper panel) and damping times 
(lower panel) of the QNMs g\,f,p\ for stellar models corresponding 
to a possible evolutive sequence of stationary configurations. 

of this point). 

Furthermore, for i — 1 the damping times of the /-mode and 
g-mode are nearly coincident, showing that in the early stages 
the g-mode is as effective as the /-mode as a source of grav- 
itational waves. Indeed, the function /gw given by Eq. d42l) . 
whose square integral over the star is the gravitational wave 
luminosity, is similar for the two modes due to a similar pro- 
file of the energy density perturbation. At later "times" x g 
becomes much larger than T/ and, as the PNS tends to the NS 
final configuration, the g-mode frequency tends to zero. 



VI. CONCLUSIONS 

In this paper we have investigated how the frequencies 
and damping times of the quasi-normal modes of a proto- 
neutron star depend on the physical quantities which char- 
acterize the stellar configurations during the quasi-stationary 
evolution. The most important is the entropy profile inside 
the star, whereas the dependence on the lepton composition is 
weaker. 

The most interesting result is that if the entropy gradient be- 
tween core and envelope is large, the frequencies of the first g- 
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mode, of the fundamental mode, and of the first p-mode tend 
to cluster in a small region near 1 kHz, whereas the damping 
time of the first g-mode and of the /-mode become compa- 
rable. This means that during the initial phases of the quasi- 
stationary evolution, when the core entropy is low and the en- 
velope entropy is large, these two modes are competitive as 
far as gravitational wave emission is concerned. 

The damping times are of the order of a few seconds, 
smaller than dissipative timescales associated to neutrino pro- 
cesses, which are of the order of 10-20 s. Thus, if the star has 
some mechanical energy to dissipate, it is likely that it will do 
it through the g\ and / modes. 
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Appendix A: The Lindblom-Detweiler equations 



The system of the LD equations ifiol [Till consists of four first-order differential equations in the quantities Hf m (r), K im (r) 
W bn (r),X hn {r): 



Hi"" = — 



£ + 1 + 



2Me . , , 



4nre A (p - e) 



K lm, = 



2r 



H, - 



£+1 w' 



K em -87t(e + p)- — W 



H^ m + K em -l67i{e + p)V e ' n 

A/2 



On 



£+1 



W em , re X/2 



X 



Cm I 



-X 



Cm 



(e+p), 



.V/2 



r 2 



(e+ P )c 2 

47t{e + p)e l / 2 + C0 2 e X / 2 -V 



(Al) 



The remaining perturbation functions, H^ m (r), V lm {r), H^ m {r), are given by the algebraic relations 



' HI I 



U-\)(t+2) o ' 

3M+- ^ '-r+Anr'p 



Cm 



£(£ + !] 



{M + Anr'p) - (0 2 r 3 e 



23 -(A +V /) 



H 



Cm 



{£ l)(*+2)_ m2 j.-y + e —(M+A%r 3 p){3M-r+An?p) 



-r— (0~r 



K 



•III 



,V/2 



X lm = co 2 (e + p)e-^ 2 V em - L/.Y-W W t«> + - — (e + p)M m , 

r 2 



Trim Trim 

H =H 2 . 



(A2) 



Equations (IAU and ( IA2b are solved numerically inside the star, assuming that the perturbation functions are non-singular near 
the center. An asymptotic expansion of the equations near r = shows that this requirement implies 



X im (0)= [e(0)+p(0)]e 



V(0)/2 



4k 

T 



[e(0) + 3 P (0)] - (o 2 —^jW em (0) + -K tm {Q) 



Hf m {0) 



£{£+l) 



2£K im (Q) + 16ff[e(0) + p(0)]W im (0) 



(A3) 



r 



On the stellar surface, r — R, one assumes continuity of pressure perturbation, i.e., 
the perturbation functions and the vanishing of the Lagrangian 



X tm (R)=0. 



(A4) 
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In the exterior, the metric perturbations are described by the 
Zerilli function 



7 tm 



J+2 



nr + 3M 



K bn -e^H 



fm \ 
1 



(A5) 



[where n = {I — 1)(£ + 2)/2], which is solution of the Zerilli 
equation 



2-7(7/ 



drl 



[co 2 -Vz(r)]Z e ">=0 



(A6) 



with = r + 2Mln(r/2M- 1) and 

_ x 2n 2 (n + l)r 3 + 6n 2 Mr 2 + \%nM 2 r + 18M 3 

V? = e t- . (A7) 

r 3 (nr + 3M) 2 V ; 



Finally, to describe free oscillations of the star we must im- 
pose the outgoing wave boundary condition 



Z bn (r)^e 



(r — )■ oo) . 



(A8) 



A solution of Eqs. (lAlb . ( IA61 > satisfying the boundary condi- 
tions ( IA3I ), (IA4I) . ( IA8b only exists for a discrete set of (com- 
plex) values of the frequency co = 2nv + i/T: the quasi- 
normal modes of the star. 
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